Production of Ultracold Molecules with Chirped Nanosecond Pulses: Evidence for 

Coherent Effects 
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We use frequency-chirped light on the nanosecond time scale to produce ultracold 87 Rb2 molecules 
in the lowest triplet state via the process of photoassociation. Comparing to quantum simulations 
of the molecular formation, we conclude that coherent stimulated emission plays an important role 
and is primarily responsible for the significant difference observed between positive and negative 
chirps. 
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I. INTRODUCTION 

Applying the concepts of coherent control to the ma- 
nipulation of ultracold systems is a topic of considerable 
current interest. Coherent control [H, [2] usually involves 
internal degrees of freedom, such as molecular vibration 
and rotation, while cooling and trapping techniques [3[ 
deal with external degrees of freedom. The time scales 
are usually quite different as well: coherent control is typ- 
ically done with ultrafast lasers while motion at ultralow 
temperatures is very slow. A particularly noteworthy 
convergence of these two fields is the formation of ultra- 
cold molecules from ultracold atoms by the process of 
photoassociation 0] (PA). This free-bound transition is 
a simple binary reaction starting with a narrow range of 
continuum energies, so coherence can be expected to play 
an important role [5j . 

In recent years, there have been many proposals for 
coherent control of PA with shaped ultrafast pulses @. 
This has been motivated in large part by the desire to 
form ultracold molecules [6( with high efficiency and state 
specificity for their many potential applications in quan- 
tum information, precision spectroscopy, ultracold chem- 
istry, and quantum dipolar systems. So far, experimental 
progress towards coherently controlled PA has been lim- 
ited to control of the photodestruction of already existing 
ultracold molecules p], H| and the observation of coher- 
ent transients in PA with femtosecond pulses [H, [l(| ■ In 
recent work, we have used frequency-chirped light on the 
nanosecond time scale to coherently control a closely re- 
lated process, laser-induced inelastic collisions. Because 
our nanosecond pulses are well matched to the long-range 
motion of the colliding atoms, we have seen that the col- 
lision rate depends not only on the chirp direction 
but also on the shape of the chirp |l2j. In the present 
work, we apply our chirped pulses to the process of PA 



and directly detect the resulting ground-state molecules. 
We see a dependence of the formation rate on chirp di- 
rection, in agreement with quantum simulations of the 
PA dynamics. These simulations reveal that despite the 
presence of spontaneous emission, a significant portion of 
the chirp dependence arises from a coherent effect: stim- 
ulated emission into a specific high vibrational level. 



II. EXPERIMENT 



In the experiment [111 ], we illuminate ultracold Rb 
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atoms with nanosecond-scale pulses of frequency-chirped 
light, forming long-range excited-state 87 Rb2 molecules 
via PA. These excited molecules subsequently radia- 
tively decay into high vibrational levels (v") of the a 
3 E+ metastable state which are detected by resonantly- 
enhanced multiphoton ionization (REMPI). The ultra- 
cold atoms are provided by a magneto-optical trap 
(MOT) operated in the phase-stable configuration and 
loaded continuously by a slow atomic beam emanating 
from a separate source MOT. The atomic temperature 
and peak density are ~150 /j>K and ~5xl0 10 cm~ 3 , re- 
spectively. 

The frequency-chirped light is produced by rapidly 
varying the injection current of an external-cavity diode 
laser with a 240 MHz arbitrary waveform generator 
(AWG) and the chirp shape is measured by a hetero- 
dyne technique. The AWG waveform is approximately a 
5 MHz triangle wave, but with programmed adjustments 
to produce approximately linear positive and negative 
chirps with equal slopes during the pulse. Each chirp 
covers approximately 1 GHz in 100 ns and is centered on 
the 87 Rb2 PA transition, a strong line located 7.79 GHz 
below the 5S 1/2 (F = 2) -> 5P 3/2 (F' = 3) asymptote 
and determined to have 0~ character [l3| . To minimize 
the amplitude modulation, the chirped light is used to 
injection lock a separate 150 mW slave diode laser [lij ]. 
A sequence of 40 ns FWHM pulses of the chirped light 
is generated by switching with an acousto-optical mod- 
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ulator (AOM). The timing selects either the positive or 
negative chirps. 

REMPI detection of the resulting Rb 2 molecules is per- 
formed with 5 ns, 4 mJ pulses from a Nd-YAG pumped 
pulsed dye laser tuned to A=601.9 nm and focused to ^3 
mm diameter at the MOT location. Based on previous 
work [TH, this light ionizes high-v" levels of the a 
3 E+ state which are expected to be populated by our PA 
to long-range excited states. The REMPI spectrum is 
similar to that from molecules produced by MOT light, 
and is dominated by a broad feature centered at 601.9 
nm. Individual high-lying v" levels are not resolved due 
to the 0.2 cm -1 pulsed laser bandwidth. The resulting 
RbJ ions are accelerated to a Channeltron ion detector 
and distinguished from Rb + ions by their time of flight. 
The timing of the experiment is as follows. A sequence 
of up to 5xl0 4 chirped (or unchirped) pulses is applied to 
the trapped atoms and 25 fis later, the REMPI pulse fires 
and the ions are detected. The entire cycle is repeated 
at 10 Hz. The MOT beams are extinguished during a 50 
/us window centered on the REMPI pulse to avoid ion- 
ization of excited atoms by the REMPI light. We use 
a train of chirped PA pulses, so in order to obtain the 
actual molecular formation rate R, we must account for 
the loss of molecules during this PA time window. There 
is slow photodestruction of a 3 E+ molecules by subse- 
quent chirped pulses at a time- averaged rate Tpp> as well 
as their escape from the detection region at a rate r esc . 
The number N of detectable molecules varies with time 
according to: 



N(t) = - ——(1-e 

1 PD + 1 esc 



We determine r esc = 108(7) s _1 by measuring the expo- 
nential decay of the REMPI signal from MOT-produced 
molecules (i.e., without chirped PA light) as the REMPI 
pulse is delayed within a fixed 7 ms window following ex- 
tinction of the MOT light. We determine Tpo and R by 
measuring the REMPI signal as a function of the length 
of the PA window (i.e., the number of chirped pulses) 
and fitting to Eq. 1 as shown in Fig. 1(a). For this case, 
Tpd ^200 s _1 for the positive chirp, implying a pho- 
todestruction probability of ^4xl0 -5 per pulse. We find 
that r p£> is linear in intensity and depends somewhat on 
chirp direction. 

The quantity of interest is the formation rate R shown 
in Fig. 1(b). These values are derived from the fits to 
curves such as Fig. 1(a) for various intensities. For each 
chirp direction, R increases with intensity, but exhibits 
some degree of saturation. The important point is that 
there is a pronounced dependence on chirp: the positive 
chirp has a rate higher than the negative chirp, but lower 
than the unchirped case. 
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FIG. 1. (Color online) (a) RbJ REMPI signals vs. photoasso- 
ciation time for unchirped, positively chirped, and negatively 
chirped pulses, along with fits to Eq. 1. The pulse repetition 
rate is 5xl0 6 pulses/s and the peak intensity of the pulses is 
fixed at 32.2 W/cm 2 . (b) Rfc>2 formation rate vs. intensity for 
the various chirps. 



III. SIMULATIONS 

To model the ultracold collisional dynamics, we solve 
the time-dependent Schrddinger equation. The dressed- 
state Hamiltonian reads: 

ff + v g j m (t) mt(t) \ 

H= hfl* (t) f + V + A , (2) 

\hni(t) o f + Vx + Aj 

where T is the kinetic energy operator and Vj (j = g, 
or 1 for the a 3 Ej, 0~ and l g electronic states) are the 
internuclear potentials computed in (l6j with coefficients 
from [l?} adjusted for the proper scattering length [l8|, 
[l9[ . The two excited states correspond to the assignment 
of 13]. A is the central detuning of the light from the 
asymptote [5S 1/2 (F = 2) + 5P 3/2 (F' = 3)]. The 0" 
detuning was shifted to yield the correct experimental 
spacing of ^0.6 GHz and to be centered on v'=78, in 
accordance with For partial waves other than s, we 
add a rotational barrier Vj = J (J + l)/2mR 2 to the a 
3 E+ potential. 

The time-dependent couplings between the a 3 E+ (g) 
and excited (j=0,l) states due to the chirped pulse are 
given by: 

mj = » gj e e [-(*=^)»4*3(t)(t-W. r )] (3) 

where fi g j are the transition dipole moments (including 
Franck-Condon factors), e is the peak electric field, a 
is the pulse width, t cen t er is the center of the intensity 
pulse, and ui(t) are instantaneous frequency offsets from 
A derived from smoothed interpolations of the hetero- 
dyne signals [20(. 



To enable efficient computation for nanosecond 
timescales, the simulations are performed within the ba- 
sis of vibrational levels calculated on a time-independent 
mapped Fourier grid [HI, [22[. A limited bandwidth is 
then taken to represent each of the vibrational Hamilto- 
nians: ~15 GHz for the 0~ and l g excited states; and 
278 GHz (16 MHz) for the a 3 S+ bound (scattering) 
manifold. We have verified that this representation is 
sufficient by extending the basis sets and checking con- 
vergence. The initial single state is a box-normalized 
scattering state at Eo = fc^T, where T = 150 /iK is the 
sample temperature. 

As described in [2Q|, spontaneous decay is taken into 
account by adding multiple sink channels corresponding 
to decay from each of the excited potentials into various 
a 3 E+ vibrational levels. Note that this model precludes 
the possibility of multiple incoherent excitations. How- 
ever, as discussed below, almost all of the population 
that decays into detectable levels is far from resonance 
and will not participate in subsequent dynamics. 

The computation gives the probability of molecular 
production, Pe ,j, for a given initial box- normalized 
state and a given partial wave, J. Representing the 
Hamiltonian (Eq. 2) on a finite grid in coordinate space 
allows us to divide the trap volume, V, into many smaller 
boxes of volume Vboxi whose number is larger then the 
number of atoms, N. Then, following [23|, the number of 
molecules, N mo i j, is computed by: 
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where n is the density, fi is the reduced mass, and ^\ E 
is the density of energy states evaluated at Eo. To find 
the molecular formation rate Pj for each peak intensity 
Io, we then multiply by the number of chirps per ms 
and spatially average over the density distribution in the 
trap and intensity profile of the photoassociation laser. 
Following (23|, we find the overall molecular formation 
rate by summing over all the partial waves necessary for 
convergence: 



P(/ ) = £(2J+l)Pj(Jo). 



(5) 



,7=0 



In Fig. 2 we plot these simulated Rb2 formation rates 
vs. peak intensity for various chirped pulses. The values 
shown correspond to the total number of molecules in a 
3 E+ (v"= 0-39) at t=200 ns after the beginning of the 
chirp. More than 93% of those molecules reside in v" =37- 
39 and thus lie within the REMPI bandwidth. The high- 
est level, v"=40, is excluded from this sum because it is 
bound by only 39 MHz and therefore easily photodisso- 
ciated by the chirped light. Also, its large outer turning 
point inhibits detection at our REMPI wavelength. The 
duration of the chirp is only 100 ns, but we allow for spon- 
taneous emission from the excited states to run its course. 
Comparing to the experimental results in Fig. 1(b), we 
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FIG. 2. (Color online) Simulated molecular formation rates 
vs. intensity for unchirped (U), positively chirped (P), and 
negatively chirped (N) pulses. The dashed curve (P') is for 
the positive chirp, but excluding the coherent contribution to 
v"=39. Also shown as the shaded region is the range of results 
for the unchirped pulses scaled by the ratio of the unchirped 
bandwidth to the chirped bandwidth (see text). The dotted 
curve U' is for a pulse shorter by a factor of 48. 



see good overall agreement, especially for the dependence 
on chirp: the rate for positive chirps exceeds that for neg- 
ative chirps, but is less than that for unchirped pulses. 

Comparing chirped and unchirped results is problem- 
atic because the pulses have different bandwidths. Ultra- 
fast pulse shaping in the frequency domain [24[ leaves the 
bandwidth fixed while stretching the pulse. In contrast, 
our addition of chirp in the time domain maintains the 
pulse width at 40 ns FWHM, but increases the bandwidth 
from the transform limit of 11 MHz FWHM to 524 MHz 
FWHM. In the simulations, we vary the center detuning 
of the unchirped pulse, with the intensity fixed at 89.29 
W/cm 2 , and find a 22 MHz FWHM in the formation 
rate. Doing the same in the experiment, we find a larger 
bandwidth of 79 MHz. The limits of the shaded region 
of Fig. 2 indicate scalings of the unchirped results by the 
ratio of each of these bandwidths to the chirped band- 
width. This scaling allows a comparison of unchirped and 
chirped results at the same intensity per unit bandwidth. 
From this point of view, both the negative and positive 
chirps are more efficient than unchirped pulses. For com- 
pleteness, we have also performed the simulations for a 
much shorter 0.84 ns FWHM unchirped pulse, increasing 
the peak intensity to keep the pulse energy fixed. The 
transform limit of this pulse is 524 MHz, the same as for 
the 40 ns chirped pulses, but its molecular formation rate 
(dotted curve in Fig. 2) is lower. This is again consistent 
with higher efficiency for chirped PA at a fixed intensity 
per unit bandwidth. 

The main conclusion to be drawn from Figs. 1(b) and 
2 is that the positive chirp gives a higher production rate 
than the negative chirp. By examining the time depen- 
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FIG. 3. (Color online) Evolution of the molecular levels dur- 
ing the positively chirped (a), and negatively chirped (b), 
pulses. The horizontal lines are the relative energies of the 
vibrational levels of the excited 0^ and l g molecular states, 
while the energies of the a 3 Ej zero-energy continuum and the 
a 3 E+ level, with the energy of the chirped photon added, are 
represented by the upper and lower black curves, respectively. 
In this picture, a curve crossing indicates when the light is 
resonant with the corresponding transition. Ground-excited 
couplings are not included in these plots. Double-ended ar- 
rows indicate the pulse widths (FWHM) . The point labeled A 
is the initial continuum state, B indicates resonance with the 
(v'=78) level, and C indicates the approach to resonance 
with the transition from 0g (v'=78) down to a 3 E^ (v"=39). 



dence of the various state populations, we have identified 
the mechanism responsible for this difference. In Fig. 
3, we show the relative energies of the excited molecu- 
lar levels involved in the chirp. These are constant in 
time. We also show the a 3 S+ zero-energy continuum 
and v"=39 level, with the photon energy added. For 
clarity, the v"=40 level, bound by only 39 MHz, is not 
shown. The time dependences of these energies reflect 
the frequency variations of the chirps: positive in (a) and 
negative in (b). At a curve crossing, the chirped light is 
resonant with the transition between the corresponding 
states. For example, at the peak of the positively-chirped 
pulse (point B in Fig. 3(a)), the chirp is resonant with 
the PA transition from the zero-energy continuum to 0~ 
(v'=78). 
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FIG. 4. (Color online) Populations of various molecular 
states during the unchirped, positively chirped, and nega- 
tively chirped pulses for I = 89.29 W/cm 2 : (a) u~ (v'=78); 
(b) l g (v'=227); (c) a 3 Ej bound levels populated by sponta- 
neous emission (SE) from 0g ; (d) a 3 E+ bound levels popu- 
lated by SE from 1 9 ; (e) a 3 E+ (v"=39) population resulting 
from stimulated emission from 0^ (v'=78). Note that in (e), 
only the positive chirp has a significant contribution. 



In Fig. 4, we plot the time-dependent populations of 
various excited and a 3 £+ states. Figs. 4(a) and 4(b) 
show the populations of the two dominant excited states, 
0g (v'=78) and l g (v'=227), respectively. As expected, 
the time ordering of population transfer to these states 
is reversed for the positive and negative chirps. The 
unchirped pulse excites only to 0~ (v'=78) since it is 
never resonant with l g (v'=227). The excited-state pop- 
ulations eventually decay due to spontaneous emission. 
As shown in Figs. 4(c) and 4(d), a small fraction of these 
decays populates a 3 E+ high-v" levels, with 0~ dominat- 
ing due to better Franck-Condon overlap. Interestingly, 
as shown in Fig. 4(e), there is another contribution to 
the v"=39 population, but only for the positive chirp. 
Referring back to Fig. 3(a), we see that towards the end 
of the positive chirp (point C), resonance between 0~ 
(v'=78) and 3 S+ (v"=39) is approached. This results 
in 0~ (v'=78) population being stimulated down to 3 E+ 
(v"=39). In contrast, for the negative chirp (Fig. 3(b)), 
this resonance is approached near the beginning of the 
chirp, when there is no excited population to be stim- 
ulated down. The time ordering of these resonances is 
crucial to the population transfer and breaks the sym- 
metry between positive and negative chirps. If we omit 
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this coherent contribution to the simulated molecular for- 
mation rate for the positive chirp, we obtain the dashed 
curve in Fig. 2, demonstrating that this contribution 
is responsible for the majority of the difference between 
positive and negative chirps. The remaining difference 
is due to the shape variation between positive and nega- 
tive chirps (Fig. 3). We have verified in the simulations 
that symmetric linear positive and negative chirps give 
the same result when this coherent contribution is omit- 
ted. We note that this coherent contribution would be 
even larger if spontaneous emission did not deplete the 
excited-state population before the stimulated emission 
occurs. 



IV. CONCLUSION 

In summary, we have investigated the formation of ul- 
tracold molecules using frequency-chirped light on the 
nanosecond timescale. We see a significant enhancement 
in the formation rate for the positive chirp relative to 
the negative chirp in both the experimental data and 
the quantum simulations. The temporal evolutions of 



the various state populations reveal the mechanism re- 
sponsible for this difference: photoassociation followed 
by stimulated emission into a high-vibrational level of 
the a 3 S+ state. Although we observed a similar trend 
(Ppos > Pneg) hi the rate constant (3 for trap- loss col- 
lisions induced by chirped light with similar parame- 
ters [ill, Ell Hil > the mechanism here is quite different. 
The collisional work utilized smaller detunings and thus 
longer-range excitation, so the time scale of the chirp and 
the atomic motion were better matched. In the present 
work, the excited state vibrational period is ~1.7 ns, 
much shorter than the chirped pulses. We expect that 
going to faster timescaies and higher intensities, together 
with controlling the details of the chirped pulses [25| , will 
allow further optimization of the molecular formation. 
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